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ABSTRACT 

To probe the behavior of polarized gluons in the proton, we propose the charmed 
hadron, such as A+ and D*, production in the forthcoming RHIC experiments. We 
found that the spin correlation between the target proton and the produced A+ 
baryon might be a good signal for testing models of the polarized gluon distribution 
in the proton. 



t Talk presented by K. Ohkuma at the XVI Particle and Nuclei International conference, 
Osaka, Sep. 30 - Oct. 4, 2002 



I. INTRODUCTION 



The advent of so-called "the proton spin puzzle" which has emerged from the 
measurement of the polarized structure function of proton g±(x) by the EMC col- 
laboration IjH, has stimulated a great theoretical and experimental activities in 
nuclear and particle physics @. Though a great deal of efforts have been made 
for solving this puzzle so far, many problems still remain to be solved. As is well 
known, the spin of proton is carried by quarks, gluons and their orbital angular 
momenta. From the next-to-leading order QCD analysis for many and precise 
data on the polarized structure function gi(x) of nucleons, now we have a rather 
good knowledge on the polarized quark distribution in the proton. However, the 
polarized gluon distribution AG(x) in the nucleon is still very uncertain. To know 
how the gluon polarizes in the nucleon is very important to solve the proton spin 
puzzle. So far, the behavior of gluons in the proton have been studied in many cases 
for deep inelastic polarized e — p scattering. However, the Relativistic Heavy Ion 
Collider(RHIC) could open another chance to probe internal structure of proton 
via pp collisions. In the RHIC experiments |J, several interesting processes, such 
as high p T prompt photon production, jet production, heavy flavor production, 
etc. are proposed to probe the polarized gluon in the proton. Here we also propose 
another processes, i.e. the polarized charmed hadron production, pp — > A+X and 
pp — > D*X, in the polarized proton-unpolarized proton collision to extract informa- 
tion about AG(x).0 In these processes, A+ is mainly produced via fragmentation 
of a charm quark which is originated dominantly from gluon-gluon fusion subpro- 
cess. Q Thus, its cross section is expected to be sensitive to the gluon distribution 
in the target proton. Moreover, since A+ is composed of a heavy charm quark and 
antisymmetrically combined light up and down quarks, the spin of A+ is basically 
carried by a charm quark which is produced via gluon-gluon fusion subprocess. 
Therefore, observation of the spin of the produced A+ gives us information about 
AG(x) in the nucleon. 

^In this report we focus only on the A+ production, though we have calculated for D* pro- 
duction, too, because the main point of the result remain unchanged. 
' 2 This is because charm quarks are tiny contents in the proton. 



II. SPIN CORRELATION ASYMMETRY AND 
ITS STATISTICAL SENSITIVITY 



As a useful observable to extract AG(x) in the proton, we introduce a spin corre- 
lation asymmetry between the target polarized-proton and produced A+ baryon; 

^ da ++ — da+- + da — da- + _ dAa/dX 

da ++ + da + - + da + de_ + da/dX ' ' 

where <icx + _, for example, denotes the spin-dependent differential cross section with 
the positive helicity of the target proton and the negative helicity of the produced 
A+ baryon. p T and 77, which are represented as X in Eq. ([!]), are transverse mo- 
mentum and pseudo-rapidity of produced A+, respectively. The spin- independent (- 
dependent) differential cross section d(A)a/dX can be calculated by the quark- 
parton model (see ref. Q for details). 

Statistical sensitivities of All for the p T and r\ distribution are estimated by 
using the following formula; 

5A LL ~ 1 1 (2) 

P V b A+ eLT ° 

To numerically estimate the value of SAll, here we use following parameters: op- 
erating time; T =100-day, the beam polarization; P =70%, a luminosity; L = 
8 x 10 31 (2 x 10 32 ) cm" 2 sec" 1 for y/s = 200 (500) GeV, the trigger efficiency; 
e = 10% for detecting produced A+ events and a branching ratio; b A + = Br(A+ — > 
pK~7i + ) ~ 5% [^]. The branching ratio of this purely charged decay mode is 
needed to measure the polarization of produced A+. o denotes the unpolarized 
cross section integrated over suitable p T or r\ region. 

III. NUMERICAL ANALYSIS 

In the numerical calculation of An, we limited the integration region of rj and p T 
of produced A+ as —1.3 < 77 < 1.3 and 3 GeV < p T < 15(40) GeV, respectively, 
for = 200(500) GeV. The range of rj and the lower limit of p T were selected 
in order to get rid of the contribution from the diffractive A+ production. As for 
the upper limit of p T , we took it as described above, for simplicity, though the 
kinematical maximum of p T of produced A+ is slightly larger than 15 GeV and 
40 Gev for ^Js =200 GeV and 500 GeV, respectively. In addition, we took the 



AAC|| and GRSV01 [0] parameterization models for the polarized gluon distri- 
bution function and the GRV98 M model for the unpolarized one. Though both 
of AAC and GRSV01 models excellently reproduce the experimental data on the 
polarized structure function of nucleons gi(x), AG(x) for those models are quite 
different. Therefore, those models should be tested in other processes. Since our 
process is semi-inclusive, the fragmentation function of a charm quark to A+ is 
necessary to carry out numerical calculations. For the unpolarized fragmentation 
function, we used Peterson fragmentation function, D c ^ A +(z) || |§. However, since 
we have no data, at present, about polarized fragmentation functions AD s ^+(x) 
for the polarized A+ production, we took the following ansatz for it 

^D^ Kt (z)=C c ^ At D c ^ At (z) } (3) 

where C C ^ A + is a scale- independent spin transfer coefficient. In this analysis, we 
studied two cases: (A) C C ^ A + = 1 (non-relativistic quark model) and (B) C C _^ A + = 
z (Jet fragmentation model [fTD| ). As we discussed before, if the spin of A+ is equal 



to the spin of the charm quark produced in the subprocess, the model (A) might 
be a reasonable scenario. 



IV. RESULTS AND DISCUSSION 

Numerical results of All are shown in Fig. [I] and Fig. 0. In these figures, we 

attached SA LL to the solid line of A LL calculated for the case of the GRSV01 

parametrization model of polarized gluon and the non-relativistic fragmentation 

model.0 Comparing with those figures, we can see that the r\ distributions of 

All are more effective than the p T distributions at y/s =200 GeV and 500 GeV 

to distinguish various models. As shown in the right panel of Fig. [| given at 

y/s =500 GeV, we could distinguish the parametrization models of polarized gluon 

as well as the models of the spin-dependent fragmentation function, though the 

magnitude of All is rather small. At y/s =200 GeV, the magnitude of All for r\ 

distribution becomes larger, though statistical sensitivities are not so small. If the 

integrated luminosity at y/s =200 GeV becomes large and the detection efficiency e 

is improved, this observable could be promising to distinguish not only the models 

of AG(x) but also the models of AD (z). On the other hand, we cannot say anything 

" 3 Note that as shown from Eq.(2), SAll does not depend on both of the model of polarized 
gluons and the model of fragmentation functions. 
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Figure 1: Spin correlation asymmetry of p T (left panel) and r\ (right panel) distri- 
bution at ^fs = 200 GeV 
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Figure 2: The same as in Fig. [T], but for y/s = 500 GeV 



from the p T distribution of A LL at high p T region, since 5A LL becomes rapidly large 
with increasing p T . However, if we confine the kinematical region in rather small 
p T range such as p T = 3 ~5(10) GeV at y/s = 200(500) GeV, it might be still 
effective. 

Though this analysis is confined to leading order, the results are interesting and 
we hope our prediction will be tested in the forthcoming RHIC experiment. 
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